UK/ 11-09 



Is 1 + Meson a Hybrid? 



Yi-Bo Yang a , Ying Chen a , Gang Li , and Keh-Fei Liu c 

Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 1000190, China 

h Department of Physics, Qufu Normal University, Qufu, 273165, P.R. China 
c Dept. of Physics and Astronomy, University of Kentucky, Lexington, KY 40506, USA 

Abstract 

_ j y 

We calculate the vacuum to meson matrix elements of the dimension-4 operator ?/>74 D iifj 
and dimension-5 operator ipEij^ji^Bk of the 1 h meson on the lattice and compare them to 
the corresponding matrix elements of the ordinary mesons to discern if it is a hybrid. For the 
charmoniums and strange quarkoniums, we find that the matrix elements of 1 ^ are comparable 
in size as compared to other known qq mesons. They are particularly similar to those of the 
2++ nieson, since their dimension-4 operators are in the same Lorentz multiplet. We conclude 
that 1 ^ is not a hybrid. As for the exotic quantum number is concerned, the non-relativistic 
reduction reveals that the leading terms in the dimension-4 and dimension-5 operators of 1 *~ 
are identical up to a proportional constant and it involves a center-of-mass momentum operator 
of the quark- ant iquark pair. This explains why 1 h is an exotic in the constituent quark model 
where the center of mass of the qq is not a dynamical degree of freedom. There are no exotics 
in QCD and hadronic models which have constituents besides the quarks to allow the qq pair 
to recoil against them. To accommodate these quantum numbers for qq mesons in QCD and 
hadron models in the non-relativistic case, the parity and total angular momentum are modified 
to P = (— ) L +'+! and J — L + I + S, where L is the orbital angular momentum of the qq pair 
in the meson. 



1 Introduction 



In the course of studying the glueball spectrum in the MIT bag model and potential 

models [jJS]) it is an underline assumption that there are valence qluons as are quarks. It is then 
a natural extension to consider hybrids of constituent quarks and gluons in the form of qqq . This 
has been studied in the potential models 0, 0], bag model fli-tlH]. flux-tube model [3, [lif . QCD 
sum rules 0, ADS/QCD [3] and lattice QCD [H-|2f3]. One of the interesting attributes of 
these hybrids is that they can have exotic J PC quantum numbers - these are J PC , s that are not 
accessible by the qq mesons in the constitute quark model where the charge and parity of a qq 
meson are given by 

P = 

C = (-) l+s , (1) 



and the angular momentum by 



J = I + S. (2) 



In light of this, these hybrids with exotic quantum numbers, particularly the 1 + has been studied 
in the above quoted references. Experimentally, there are two candidates for the 1 *~ - one is 



7Ti(1400) [27| and the other is 7Ti(1600) [28|]. They are observed in the rjn and pn channels. 

In view of fact that exotic quantum numbers are not accessible by the constituent quark- 
antiquark pair, it is suggested that the interpolation field for the hybrids of the qqg type will 
necessarily involve a gauge field tensor, i.e. of the form ipQipG where O involves 7 matrices and 
covariant derivatives and G stands for the field tensor G^ u . It is an operator with dimension > 5. 
However, it was pointed out by B.A. Li more than 30 years ago that these exotic quantum number 
can be constructed from the quark bilinears ^Qijj without the field tensor [29(. For example, 
the J PC of ^74 D iip is 1 ^ which is a dimension-4 operator. This type of operators have been 



constructed on the lattice 22| and lattice calculations have been calculated with them instead of 
the dimension-5 operator Eijk'ipjjipBk 22^, 2^- The exotic mesons can be in the form of tetraquark 
mesoniums qqqq which will require a dimension 6 interpolation field. We will not address them in 
the present work. 

The existence of the dimension-4 operator for 1 ^ that does not involve the gauge filed tensor 
raises several questions: 

• Since there exists an interpolation operator which does not involve the field tensor, does that 
mean the meson with this interpolation field is not an hybrid? One could point out that the 
dimension-4 operators involve a covariant derivative which allows it to couple to a constituent 
gluon, unlike the dimension-3 operators ^Ttp, where T is a 7 matrix, for the pseudoscalar and 
vector mesons. However, one can counter this argument by pointing out that the tensor meson 
(2 ++ ), like 1 , does not have dimension-3 interpolation field. The minimum dimension of 
its interpolation field is a dimension-4 operator ipji D jtjj which is very similar to that of the 
dimension-4 operator for 1 ^ [29[] and yet 2 ++ is an ordinary quantum number. So, how does 
one find out whether a meson is a hydrid or not? 

• Since 1 ^ is an exotic quantum number, how come one can have an operator which does not 
involve the field tensor? If one carries out a non-relativistic reduction of the operator, would 
one be able to reveal why it is not accessible to the constituent quark model? 
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To answer these questions, we shall establish criteria for identifying the hybrid and carry out 
a lattice calculation with both the dimension-4 and dimension-5 interpolation fields to compare 
their respective spectral weights against those of ordinary mesons. We will also carry out a non- 
relativistic reduction to figure out why the exotic quantum numbers are not accessible to the 
constituent quark model. We shall present the meson interpolation fields organized in dimensions 
3, 4, and 5 for various mesons in Sec. II, set criteria for distinguishing hybrids from ordinary 
mesons, and discuss the origin of the exotic quantum numbers. The numerical details are given in 
Sec. Ill and the results are given in Sec. IV. We will end with a summary in Sec. V. 



2 Formalism 

We shall discuss several types of meson interpolation fields and set up criteria in order to distin- 
guish the hybrids from the ordinary mesons via the vacuum-to- meson transition matrix elements. 



2.1 Meson interpolations fields and criteria for hybrids 

In lattice calculations, one relies on interpolation fields with the desired quantum numbers 
(e.g. J PC , isospin, strangeness, etc.) to project to the physical spectrum with the corresponding 
quantum numbers. In the following we give a list of these interpolation fields for the low-lying 
ordinary mesons (pseudoscalar, vector, axial- vector, scalar and tensor) and 1 '". They are classified 
according to the following types: 



ipTip (r is a gamma matrix), a dimension-3 operator, is labeled as the T-type; 
ipT x Dip [D = D — D), a dimension-4 operator, is labeled as the D-type; 
ipF x Bip (Bi = T;£ijkGjk), a dimension-5 operator, is labeled as the B-type. 



A more complete list can be found in Ref. 22]. 



Table 1: Interpolation operators ipYip (dimension 3, T-type), ipT x Dip (dimension 4, D-type), and 
x Bip (dimension 5, B-type). £j = ^SijkO'jk and repeated indices are summed over. 
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Here, we only list 1 + as an example of mesons with exotic quantum numbers that cannot 
be accessed by dimension-3 operators. We should point out that ordinary J = 2 mesons do not 
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have dimension-3 interpolation operators either. There are two kinds of dimension-4 1 + operators 
("0 a 74 -D itp a and e«^'0 a Sj D k4> a )- These two kinds of operators have very similar non-relativistic 
forms as will be discussed in Sec. 



A meson correlator at zero momentum is 

C ij (t) = Y / (O i (x,t)O J (0,0)). (3) 

X 

At large time separation, it is dominated by the lowest state of the spectrum with the prescribed 
quantum number 

Cij(t) __> -L(0|O i |M)(M|O J -|0)e-- i (4) 

where m is the mass of the lowest state. Besides the mass, one also obtains the vacuum to meson 
transition matrix elements (0|O|M). 

We should point out that, notwithstanding claims in many lattice calculations, the interpolation 
operators do not necessarily reflect the nature of the composition of the hadrons. They merely reveal 
how strongly the operators couple to the specific hadron, such as realized in decay constants. For 
example, the topological charge operator G^G^ projects to r\ and r/ strongly. From the anomalous 
Ward identity for massless fermions d^A® = -^jG^G^u, one has 



(0|^<V<W> (5) 



This does not mean that r/ is a glueball, even though the matrix element is larger than the matrix 
element of the isovector axial- vector current for the pion 

(0|a M Aj|7r)=m2/ w , (6) 

due to the larger rj' mass as compared to pion. In fact, the flavor-mixing angle between 771 and rjs 
for r/, r/ have been well studied with the help of axial anomaly [3(j. Including the glueball mixing 
from the KLOE experiment of <p — > JVi iv' '> the matrix elements of (0| j^tG ^ V G ^ v \M) for M = rj,rj' 
and glueball G are found to be of the same order, even though in the large N c analysis, the matrix 



elements for 77,77' are parametrically smaller by 0(l/\/N c ) than that of the glueball [311 ] . This is 
known to be related to anomaly. On the other hand, the matrix element (0| 975?! G) is more than an 
order of magnitude smaller than those of {^\q^q\r}, r]')(q = u, d, s) [31]. This shows that the lower- 
dimension quark field operators couple to the glueball much weaker than to the qq mesons. This has 
been taken as a criterion to distinguish the glueball from the qq mesons under the condition that 
the glueball does not mix with the qq mesons strongly. In view of this analysis of the pseudoscalar 
mesons, it is suggested (32| that the smallness of the matrix element of lower-dimension quark 
operator compared to those of established qq mesons is a better signal for the glueball than those 
with the higher dimensional glue operators. By the same token, we shall adopt a similar criterion 
for detecting the hybrids by examining the dimension-4 D-type matrix element (0|t/T x D ip\M) 
and the dimension-5 B-type matrix element (0|t/T x Bip\M) of the 1 ^ and compare them with 
those of the other ordinary mesons. If the D matrix element of 1 h is much smaller than others 
and the B matrix element much larger than (or at least as large as) the others, then it is a hybrid. 
Otherwise, it is not. Special attention will be paid to the comparison with the 2 ++ meson. Neither 
1 ^ nor 2 ++ has dimension-3 interpolation field and their dimension-4 operators are in the same 
Lorentz multiplet i.e. D u ip. 



2.2 Non-relativistic Operators 



To address the question of the exotic quantum number, it would be useful to find out the non- 
relativistic form of the interpolation operators listed in Sec. 12.11 We use Foldy-Wouthuysen-Tani 
transformation [i^] for non-relativistic reduction to the heavy quark and anti-quark fields described 
by the Pauli spinors <ft and x- 



Table 2: Non-relativistic form for the three kinds of operators (T, D and B) as shown in Table 1. 
Here we list the operators O in the interpolation field x^^cf). Repeated indices are summed over. 
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The Dirac spinor ^ and ^ are expanded in terms of cj> and x in 1/m as 
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(7) 
(8) 
(9) 

(10) 



Operator D used here is the spatial part of the covariant derivative and m is the heavy-quark 
mass. The Pauli spinors 4>/4>^ and x /x are the annihilation/creation operators for the heavy quark 
and antiquark which satisfy the relation 



0|O) =X f |0) =0; (O|0t = ( | x = 0. 



(11) 



With the above approximation, we could reduce the operators listed in Table Q] with a given 
J PC to the form of an d f O'X with O now involves cr, D , and B. We shall still classify 

them according to their dimensions and label them the same as before, i.e. T-type (dimension 3), 
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D-type (dimension 4), and i?-typ e ^d imension 5). The operators for y^O(\) to leading order in 1/m 
are listed in Table [2J Note that D acts on the quark and anti-quark fields, while d acts on the 
glue field B. 

2.3 Exotic quantum numbers 

We see from the non-relativistic reduction in the above section that the dimension-4 (.D-type) 
interpolation field for the 1 ^ meson involves a symmetric combination of *5 and ll This is the 
center of mass momentum operator of the qq pair. We now see why this operator is not admissible 
in the quark model with only the constituent quark degree of freedom. In this model, the center of 
mass of qq is not a dynamical variable due to translational invariance, while the quantum number 
J PC is defined in the center of mass of the qq pair. In QCD, there are gluons besides the quarks 
so that the qq pair can have orbital angular momentum relative to the glue stuff, much like the 
orbital motion of the electron pairs around the nucleus in the atom, or the planetary motion of the 
earth-moon pair around the sun. This is also true in models where there are other constituents 
that the qq pair can recoil against. For example, in the MIT bag model, the qq can have orbital 



angular momentum against the bag if the latter is made dynamical 36] . In the chiral quark model, 
the qq can recoil against the pion and so can it against the flux tube in the flux-tube model. 

In fact, the issue of the the exotic quantum number and its relation to the center-of-mass motion 
of the qq has been raised in the MIT bag model [l|, [35j. An example is given for the 2 +± meson 
where the quark and anti-quark orbital wavefunctions are given as 

Since both C = ± are possible, they double the spectrum from the conventional constituent qq 
model. It is pointed out that the symmetric combination leads to a P-wave for the center of 
mass of the qq. In the nuclear shell model with harmonic oscillator potential, this is considered a 
spurious center-of-mass excitation since the center of mass is pinned down by the harmonic oscillator 
potential. If the bag is not dynamical like the external harmonic oscillator in the shell model, it 
can be removed with center-of-mass correction [37}]. However, if the bag is considered dynamical 
with surface fluctuations [3(1], this center-of-mass motion is physical and so is the 2" 1 state. By 
the same token, one can consider the 1 _± states with the combination 

*(1 _± ) = -=(1Si2Si T2SiLSi), (13) 
^/2 V 2 a 2 a" 

with the anti-symmetric combination being the 1 ^ state where both the center-of-mass and relative 
coordinates are in the P-wave for harmonic oscillator wavefunctions. 

To conclude this part of the discussion, we see that there are no exotic quantum numbers in 
QCD and models with additional constituents besides the qq pair. The 'exoticness' is only in the 
context of the constituent quark model with only qq degree of freedom. These quantum numbers 
can be accommodated with parity and total angular momentum from Eqs. (JTJ) and ([2]) supplanted 
by 

p = (_) L + i+1 

J = L + l + S. (14) 
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where L is the orbital angular momentum of the qq pair in the hadron. The charge parity 
C = {—) l+s remains the same, provided that other degrees of freedom in the hadron are not excited 
and gives C = +. In the case ofl ^,L = l = S = l and the two operators in Table [2] correspond 
to S + L = and S + / = respectively. Other exotic quantum numbers, e.g. 0" 1 , 2 H ,3 h can 
all be accommodated in Eq. (Q2 



3 Numerical Details 

We shall give lattice details including the action, the parameters as well as the operators used 
for the interpolation fields of various mesons. 

3.1 Improved Clover Action 



We adopt the anisotropic Wilson gauge action 38] in the quenched approximation for the 



A X y &xy 



present study. The improved anisotropic Wilson fermion action is 

M-xy — $xy® ~t~ Axy 

r 3 

1/ (2/Wr) + Pt ^2 COi-^Oi + Pa{p\lT\1 + 023-^23 + C31.F31) 
i=l 

~ ^2vn[{l-l f i)U f i(x)5x + L l ,y + (l + "f fl )Ul(x- [i)5 x -n,y , ( 15 ) 
n 

where the coefficients are given by 

v £11 



Vi = 7T- ' % = o > a 



2u ' ' 2 7 2 



max 



(1 + v 

Pt = . p- = • (16) 

with ^ = a s /a^ being the bare aspect ratio of the asymmetric lattice, and v the bare speed of light 
parameter. Another parameter u s , taken to be the fourth root of the average spatial plaquette 
value, is used to incorporate the tadpole improvement of the spatial gauge link Ui(x). 
With this fermion action, the bare mass of the quark is 

m oa s = x £ - 3^ . (17) 

The lattice used in this study is of the size 12 x 12 x 12 x 96 at (3 = 2.8 which gives a s = 0.0276 
fm, with the aspect ratio £ = a s /at = 5. 

The bare k of the charm quark is set to 0.060325 with the bare speed of light parameter v = 0.74, 
which is determined by fitting the mass of J/ip. Similarly, The bare n of the strange quark is set 
to 0.0615 which gives the vector mass close to that of <p. 



3.2 Masses and Vacuum to Meson Transition Matrix Elements 

To construct two-point functions, we use the T-type wall operators for mesons which have 
dimension-3 interpolation fields. For those which do not have dimension-3 interpolation fields, we 
use the B-type wall source to enhance the signals. This is illustrated in Fig. CO B w denotes the wall 
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f t=o t t=0 




Figure 1: Sketch of two point functions. The lines denote quark and antiquark propagators. The 
black dot is the glue field tensor B attached at the quark wall source and sink. 



source and sink for the -B-type operator. We note that the glue field tensor B can be attached to 
either the quark field or the antiquark field at the source and sink. It is indicated by a black dot in 
the figure. When both the wall source and sink are of the -B-type operators, it is necessary to sum 
the two kinds of diagrams (middle one in Fig. [1]) to obtain an eigenstate of charge parity. When 
the sink is the point operator which has a definite charge-parity, one diagram with the B attached 
to either the quark or anti-quark wall source (the right diagram in Fig.. [T]) will suffice. T p , D p and 
Bp denote the point sinks. 

The wall-source is placed on 16 of the total 96 time slices separately for each of the 1000 
configurations to gain statistics. We calculate the correlators with both the source and sink being 
the wall B-type operators (B w ) and with -B-wall and point sinks with the T,D and B operators 
(T„,D P and B„). The color magnetic field B is smeared twice for the wall source and sink and the 
double antisymmetric derivative operator Eij^tp D j D ^ip is used to replace ifiipBi for the point sink. 

The ground state mass and the vacuum to meson transition matrix element are extracted from 
the following correlators: 

<0(T >< mLl(t) 0(T X fl^UKO)) >N V \ { ° mT X y WI^ PC )| 2 {e -mt + e -m(n T - t)] 



W) x B )^,(°))--- W 0|O ' |jPC>(JfC f (r x m ^ 0) (e- + 

(18) 



where Ny = L 3 is the three-volume factor. From these two equations, one can obtain the matrix 
element (0\O p \J PC ). 

(0\O P \J PC ) = I [ J > f 2 (O p mr x B^Uu). (19) 

N v ([if)(T x B)^)y wall [V?(r x B)i;} wall ) 

Similarly, we also obtain the masses and T-type matrix elements with the T-type wall source (right 
and left diagrams in Fig. [1]) . 
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4 Numerical results and Discussion 



4.1 Charmoniums 

We first calculate the masses and the matrix elements for the charmonium with the charm quark 
k = 0.060325 that was tuned to the physical J/^f mass. The masses from different correlators are 
listed in Table O 

Table 3: Masses of charmonium states from T- and -B-type sources and point sinks. 
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The effective masses of rj c (0 and J/^(l ) are plotted in Fig. [2] with the B- and T-type wall 
sources and T p ,D p and B p for the zero momentum point sinks. The effective masses of Xco(0 ++ ) 
and Xci(l ++ ) from the wall sources are plotted in Fig. [3] for several point sinks. The effective 
masses of Xc2(2 ++ ) and r/ c i(l '") from the wall sources are plotted in Fig. [5]for several point sinks. 
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Figure 2: Effective mass plot for rj c and J/vP with B- and T-type operators as the wall sources and T p , D p 
and Bp for the zero momentum point sinks. 

We see that the masses obtained from different correlators with different sources and sinks 
are all consistent with each other and the pattern of the charmonium masses, besides 1 are in 
reasonable agreement with experiments, except the Jryperfine splitting which is known to be smaller 
than experiment for the quenched approximation [40|. We note that dimension-4 and -5 operators 
produce the same mass of 1 ^ within errors. We take this to imply that they are the same state. 

The effective masses of the pseudoscalar (rj c ) and vector (J/*) charmonium are plotted in Fig. [2] 
for the cases with T w and B w sources and T p , D p and B p sinks. The effective masses for the scalar 
(Xco) and axial- vector (xci) are plotted in Fig. [3] and those for the tensor (xc2) and 1 h are plotted 
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Figure 3: The same as Fig. [2] for Xco and Xd 
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Figure 4: The same as Fig. [2]for Xc2(2 ++ ) and fy c i(l + ) 



in Fig. [U As we can see from Table they agree for different sources and sinks within errors. 
Different interpolation fields project to the same lowest states in all channels studied here. 

Before we discuss the results on the matrix elements, we should point out a relation between 
the dimension-4 L>-type and dimension-5 B-type operators in the non-relativistic limit. 

The double derivative operator in the leading non-relativistic expansion of the 1 ^ interpolation 
field can be expanded as 

= xH®i'dj ~ ~ di.dj + dj.di)il) 

+X i (2i.d j (A i ) - 2[A i ,A j ])il> + xH2ift) i .A j + %.A^ 



(20) 



Since we are projecting to the zero momentum meson state in the lattice calculation with periodic 
condition in the spatial direction, we have 



J d 3 x x^*d i-*d jip = J d 3 xx^~$ J d z x\^ i- d jip = J d 3 xx**dj-~$iip, 



d 3 xdi( X ^M) = I d 3 x X H i d l -A j + Aj.HM + x^dj(Ai)i> = 0. 



(21) 
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From Eqs. (JH and (j2lj) . we obtain 

J d 3 xxH^i^j + ^j~i>i)i> = J d 3 x X \2i.d 3 (Ai) - 2i.di(Aj) - 2[A i ,A j ])i/> 

d 3 x 2ix^Gijip (22) 
Therefore the leading non-relativistic terms of the zero-momentum 1 ^ interpolation fields are 

y at.k. 2m y 

2m/ dZx ^2a j G ji (t) = -^2e ijh j d 3 x^(J 3 B k <\) (23) 

fd 3 x$ a £ ijk X j < L) k i> a , ^ [ d 3 x X HDi(T- i D l5i)<j> 

J N.R. 2m J 

= 7^ J d 3 xx^2a J G ij <P = ^2e ijk J d 3 x X ^ J B k ^. (24) 

We see that, up to a sign and a proportional constant (i.e. heavy quark mass m), both 
dimension-4 D-type operators of 1 ^ are equivalent to the dimension-5 -B-type operator with 
the magnetic field in the non-relativistic limit. The matrix elements from all three operators are 
expected to be the same up to a known constant and 0(a) for heavy quarkoniums. 

The matrix elements for the charmoniums are listed in Tabled! 



Table 4: Matrix elements < OlCUJ > for charmoniums. 
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0.044 ±0.003 


0.00080 ± 0.00008 


1- 


-+ 




0.0059 ± 0.0005 
0.0054 ± 0.0004 


0.0082 ± 0.0006 



The matrix elements of 0.0059(5) and 0.0054(4) for the two D operators of 1 h are the same 
which are expected from the above discussion. Eq. ()23[) also shows that they should be 1/ma times 
that of the dimension-5 U-type to O(o). On the anisotropic lattice used here, dimension-less 
should be replaced by the anisotropic form, 

— =>- ~ 0.7048 (25) 

ma m c atq 

where £ and v are defined in Eq. (|15p . Multiplying this factor to the £>-type matrix element 
0.0082(6) gives 0.0058(4) which agrees with the -D-type matrix elements quite well. 

Furthermore, comparing T and D operators for the P-wave states ++ , 1 ++ and l" 1 in Table [2] 
shows that they are related by Thus, we expect dimension-3 T matrix elements to be e = t^t = 

10 



0.3524 times the dimension-4 D matrix elements. To check this, we plot e times the T m.e. against 
the D matrix elements in Fig. 14. II for these states and also 2e times the D matrix elements of 1 *~ 
meson against the corresponding B matrix elements. We fit the ratio of all the data and find the 
slope to be 0.35(4). This is quite consistent with e = 0.3524. This shows that the matrix elements 
we studied for the charmonium states are quite non-relativistic in the sense that higher orders in 
1/m are not important to spoil the equivalence relation we found in Eq. (j23|) and that cutoff effect 
in O(o) is small. There is a concern about mixing between operators with different dimensions. 
The results in Fig. 14.11 suggest that the mixing effects are also small. 
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Figure 5: Global fit for the ratios of T for Xco, Xa and h c (2 D for 1 + ) m.e. to the corresponding m.e. of 
D (B for 1~+). 

At first sight, the D matrix elements of 1 h are about an order of magnitude smaller than 
those of the other mesons. However, upon comparing with 2 ++ in Table we see that the 1 
operators have an extra factor of ()5 + ~D)/2m which is the velocity of the cc pair. Since the speed 
of the charm quark in J/ty is about 0.3 c, we estimate the extra factor to be ~ 0.3 (and likely to be 
less). Dividing this factor from the 1 ^ D matrix elements gives ~ 0.20(2) which is about a factor 
of two from that of the 2 ++ meson and comparable in size to the matrix elements of the other 
charmonium states. Since the matrix elements of the lowest dimension operators (i.e. D-type) of 
the 1 ^ in the charm region are comparable to and mostly smaller than those of the other known 
charmonium states, it is not a hybrid by the criteria discussed in Sec. 12.11 On the other hand, the 
matrix element of B for 1 ^ is comparable to those of the other charmoniums, except 2 ++ which 
is an order of magnitude smaller. This is presumably due to the factor of D /2m in the 2 ++ B 
operator in Table [2j Incorporating this factor of ~ 0.3 brings B matrix element of 2 ++ to within 
a factor of 3 from that of the 1 K In fact, all the P-wave operators have this D /2m factor and 
their matrix elements will be comparable or larger than that of 1 ^ when this factor is taken into 
account. The fact that 1 ^ does not have an extraordinarily large B matrix element compared to 
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other known charmonium states enhances the notion that it is not a hybrid in the charm region. 

4.2 Strange quark mesons 

Next, we consider lighter quarkonium with the strange quark. The strange meson (ss) masses 
in MeV are listed in Table El 

Table 5: Masses of strange quarkoniums from r^- and -B^-type sources and point sinks. 







B w y Tp 


B w —?■ B p 


B w — > Dp 


o-+ 


714 ± 9 


750 ± 15 


713 ± 9 


714 ± 10 


1~ 


1027 ± 9 


1030 ± 12 


1030 ± 15 


1024 ± 12 


0++ 


1570 ± 63 


1566 ± 21 


1568 ± 21 


1567 ±21 


1++ 


1580 ± 35 


1562 ± 21 


1597 ± 40 


1522 ± 39 


1+- 


1613 ± 35 


1569 ± 18 


1608 ± 54 


1598 ± 19 


2++ 






1638 ± 21 


1611 ±60 


1-+ 






2066 ± 62 


2115 ±85 



The effective masses of ?7 S (0 
plotted in Fig. El El and El 



:i""),/o( s )(0 ++ ),/i (s )(l ++ ),/ 2 ( s )(2 ++ ) and the ss 1~+ are 
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Figure 6: Effective mass plot for rj s and 4> with B w - and T^-type wall sources and T P ,D P and B p for the 
zero momentum point sinks. 



We see from Tableland Figs. [6l O and [8] that the masses from different sources and sinks are 
the same within errors. The matrix elements for < 0\O p \ J PC > for the ss mesons are listed in 
Table El 

For the light quarkonium ss, we do not expect the non-relativistic equivalence between the D- 
type and -B-type operators to hold. We shall compare the matrix elements directly. It is worthwhile 
noting that the dimension-4 D matrix elements of 1 is comparable to that of the 2 ++ meson and 
are not particularly smaller than those of the other ss mesons. Although the B matrix element of 
1 ^ is larger than that of 2 ++ , but it is not larger than those of other mesons. From these data, 
we see no evidence to distinguish the 1 ^ ss from other established ss mesons and identify it as a 
hybrid. 
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Figure 7: The same as Fig. [5] for f Q ^ and fi 
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Figure 8: The same as Fig. [6] for / 2 ( s ) and ssl 



5 Conclusion 

_ ^ y 

We set out to address the question: in view of the fact that there is a dimension-4 ^74 D ip 
interpolation field for 1 which does not involve the gauge field tensor, how does one identify it 
as a hybrid and distinguish it from the ordinary mesons, which also have dimension-4 interpolation 
fields with a covariant derivative and dimension-5 interpolation fields involving explicitly the color 
magnetic field B in the form of eijk'ip'jj X Bkip! We emphasize that one cannot judge the nature 
of a state by the appearance of its interpolation field. This is amply illustrated by the large 
matrix element (0\GG\r], rf) which shows that even though ij and rf can be produced with the glue 
interpolation field, it does not mean that they are glueballs. The glueball nature will be better 
revealed by a weak coupling to the qq interpolation field. We have also come up with an example 
where the zero momentum operators of eijk^lj x B^ip and ^74 D ip for the heavy quarks are the 
same up to a proportional constant, which is the quark mass. This implies that the former operator 
with a field tensor does not necessarily project to an excited glue state, it could project to a state 
with the qq pair in a P-wave in the hadron as the latter interpolation field in the non-relativistic 
limit suggests. 

In light of this, we compare the matrix element of (0^74 D tp\l '") and {0\ip£ijk'yjipBk\l to 
the corresponding matrix elements of the other known qq mesons. In the case of charmoniums, we 
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Table 6: The matrix elements < 0\O p \J PC > for strange quarkoniums. 





r p 


D 


B 


o-+ 


0.0247 ± 0.0002 


0.021 ±0.002 


0.005 ± 0.0001 


1 — 


0.0141 ± 0.0002 


0.0113 ±0.0005 


0.0025 ± 0.0001 


0++ 


0.043 ±0.006 


0.033 ±0.005 


0.017 ±0.004 


1++ 


0.029 ±0.004 


0.034 ±0.004 


0.0018 ± 0.0002 


1+- 


0.019 ±0.006 


0.029 ±0.005 


0.0019 ± 0.0004 


2++ 




0.010 ±0.007 


0.0003 ± 0.0001 


1-+ 




0.007 ±0.001 


0.004 ±0.001 






0.006 ±0.002 





find both the D- and -B-type matrix elements of 1 ^ are about the same size as the other mesons. 
When a velocity of the cc pair is taken into account, they are also comparable to those of Xc2(2 ++ ), 
which is most similar to 1 h in that neither has dimension-3 operator and their dimension-4 
operators are in the same Lorentz multiplet. We have also examined the strange quarkoniums and 
found that the D- and -B-type matrix elements of 1 h are comparable in size to those of the other 
ss mesons. Based on these data, we conclude that there is not much distinction between 1~ + and 
other known qq mesons. There is no evidence for it to be a hybrid. 

The leading non-relativistic expansion reveals that the dimension-4 operator 1 ^ takes the form 
of X^2m _ ( CJ ' D D i + D i<j- D)(t> and 2^r 0^ i a j D j + aj D jDi)(f>. They involve a P-wave of the qq 
pair. Since the center of mass of the qq in a constituent quark model is only a kinematical degree 
of freedom, confined center- of-mass motion is not admissible in the constituent quark model. This 
is why the J PC of 1 h and others involving the angular momentum of the qq pair are considered 
'exotic'. 

In QCD, the qq pair can recoil against the glue in the meson. Similarly, qq pair can have orbital 
angular momentum relative to the bag in the MIT bag model, to the pion in the chiral quark 
model and to the flux-tube in the flux-tube model. Thus in QCD and in models with additional 
constituents other than the qq pair, there are no exotic quantum numbers. These additional 
J PC quantum numbers can be accommodated by supplanting the parity and angular momentum 
relations to P = (-) L+l+1 and J = L + 1 + S. 
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